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We investigate the quantum Zeno effect in the case of electron tunneling out of a quantum 
dot in the presence of continuous monitoring by a detector. It is shown that the Schrodinger 
equation for the whole system can be reduced to Bloch-type rate equations describing the combined 
time-development of the detector and the measured system. Using these equations we find that 
continuous measurement of the unstable system does not affect its exponential decay to a reservoir 
with a constant density of states. The width of the energy distribution of the tunneling electron, 
however, is not equal to the inverse life-time - it increases due to the decoherence generated by 
the detector. We extend the analysis to the case of a reservoir described by an energy dependent 
density of states, and we show that continuous measurement of such quantum systems affects both 
the exponential decay rate and the energy distribution. The decay does not always slow down, but 
might be accelerated. The energy distribution of the tunneling electron may reveal the lines invisible 
before the measurement. 
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I. INTRODUCTION 

The interplay between quantum dynamics and quan- 
tum measurements has attracted much attention of 
physicists since the discovery of quantum mechanics. In 
fact, it was suggested that frequent or continuous obser- 
vations of an unstable quantum system can inhibit or 
slow down its decay This phenomenon is known as 
the quantum Zeno effect. Usually this effect is associ- 
ated to von Neumann's postulate in the theory of quan- 
tum measurements ||^. Indeed, in the classical example 
of two-level systems, the probability of a quantum tran- 
sition from an initially occupied unstable state is given 
at short time, Ai, by -P(At) = a(Ai)^. If we assume 
that Ai is the measurement time, which consists in pro- 
jecting the system onto the initial state, then after N 
successive measurements the probability of finding the 
unstable system in its initial state, at time t = NAt, 
is Q{t) = [1 - a(At)2](*/'^*). It follows from this result 
that Q{t) — > 1 for At — > 0, i.e. suppression of quantum 
transition. 

During last years the Zeno effect has become a topic 
of great interest. It has been discussed in the areas of 
radioactive decay polarized light [Q, the physics of 
atoms , neutron physics 0| , quantum optics |^,|| and 
mesoscopic physics [p^-p^ . As a matter of fact, the theo- 
retical and experimental efforts has been mainly concen- 
trated on quantum transitions between isolated levels p^ ] 
characterized by an oscillatory behavior. In this latter 
case the slowing down of the transition rate has, indeed, 
been found. However, this was attributed to the deco- 
herence generated by the detector without an explicit 



involvement of the projection postulate ]l^ , p^ . On the 
other hand, the slowing down of the exponential decay 
rate still remains a controversial issue, despite the fact 
that it is extensively studied p7|-^ and further investi- 
gations are clearly desirable. 

In this paper we focus our attention on observation of 
spontaneous decay using a microscopic description of the 
measurement device (detector). The latter point should 
be very essential in any investigations of measurement 
problems. The quantum-mechanical description of the 
measurement device would allow us to study thoroughly 
the measurement process without explicit use of the pro- 
jection postulate, or introducing different phenomenolog- 
ical terms in quantum evolution of the measured system. 
Yet, the detector is a macroscopic system, the quantum 
mechanical analysis of which is very complicated. Thus 
one can expect that mesoscopic systems, which are be- 
tween the microscopic and macroscopic scales, would be 
very useful for this type of investigation In addi- 

tion, the actual experimental investigations of the quan- 
tum Zeno effect in mesoscopic systems are within reach 
of nowadays experimental techniques. 

Our quantum-mechanical treatment of the entire sys- 
tem, including the detector is based on the new Bloch- 
type rate equations, which we derived here from the 
many-body Schrodinger equation. These equations 
would allow us to trace the quantum-mechanical behavior 
of the entire measured system during the measurement. 
In contrast with the standard Bloch equations, our equa- 
tions describe the reservoir states too. As a result, we 
can find the influence of the measurement on the energy 
distribution of the decayed system, as well as its time- 
evolution. 
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The paper is organized as follows: In Sect. 2 we gave 
a general quantum-mechanical description of an unstable 
system. We concentrated on conditions under which an 
exponential and non-exponential decay can be obtained, 
yet without including any measurement apparatus. The 
latter is introduced in Sect. 3 where we derive the rate 
equations, describing the dynamical evolution of the mea- 
surement process. The results of our analysis are de- 
scribed in Sect. 4. The last section is a summary. Some 
of the results described in this paper were presented in 



II. DYNAMICS OF A DECAYED SYSTEM 

Let us consider a mesoscopic quantum dot coupled to 
an electron reservoir on its right. Fig. 1. We assume only 
one level, Eq, inside the quantum dot, Fig. 1. The sys- 
tem is described by the following tunneling Hamiltonian 
written in the occupation number representation: 



Hqd — -EoCqCq 



E 



H.c.]. (1) 



Here the operators cJ(cq) correspond to the creation (an- 
nihilation) of an electron in state a in the reservoir, and 
the operator cj(co) is similarly defined for the state in 
the quantum dot. The fla are the hopping amplitudes 
between the states Eq, Ea. (We choose the gauge, where 
JIq. are real). This Hamiltonian (|^) does not include any 
coupling with radiation or phonon fields as had been done 
for example in 
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Fig. 1: The energy level Eq of the dot occupied by an 
electron, which tunnels to continuous states Ea of the 
reservoir. 

The initial state of the entire system cj|0 > corre- 
sponds to occupied quantum dot and empty states in 
the reservoir. This state is unstable: the Hamiltonian 
(H) requires it to decay to the continuum states cJ^jO > 
having the electron in the reservoir. The electron wave 
function, describing its evolution can be written in the 
most general way as 



6o(t)4 



{t)c 



|0) 



(2) 



where bi(t) is the time-dependent probability amplitude 
to find the system in the corresponding state i with the 
initial condition 6o(0) = 1 and ba(0) = 0. Substitution 
of Eq. (^) into the time-dependent Schrodinger equation 
idt\'^{t)) — H\'^{t)), leads to an infinite set of coupled 
linear differential equations for the amplitudes b{t). Ap- 
plying the Laplace transform 



b{E) = / b{t) exp{iEt)dt 



(3) 



and taking into account the initial condition, we trans- 
form the differential equations for b{t) into an infinite set 
of algebraic equations for the amplitudes b{E): 



(E - Eo)boiE) -Y,^o.ba = I , 

a 

{E - Ea)baiE) ~ QaboiE) = 



(4a) 
(4b) 



In order to solve these equations we replace the amplitude 
ba in Eq.(^a|) by its expression obtained from Eq.(^). 
Then we obtain 



E-Eo~Yl 



E — Ea 



ME) = i. 



(5) 



Since the states in the reservoir are very dense, one can 
replace the sum over a by an integral over Ea- 



„ E — Ea 



n^{Ea)piEa 
E — Ea 



-dEa 



(6) 



where p{Ea) is the density of states in the reservoir. To 
evaluate this integral, we can split the integral into its 
principal and singular parts, —iS{E — Ea)- (Notice that 
E = E + ie in the Laplace transform, Eq. (^). As a 
result the original Schrodinger equation (0) is reduced to 



E- Eq- A{E) + i 



boiE)=i, 



{E - Ea) ba{E) - n{Ea)bo{E) = 0, 



(7a) 
(7b) 



where r{E) = 2T:p{E)nl{E) and A{E) is the energy- 
shift due to the principal part of the integral. 

Using Eqs. (j^) and the inverse Laplace transform one 
can obtain the occupation probabilities of the levels Eq 
and Ea Yet, Eqs. (^ are not convenient if we wish 
to include the effects of the environment. These effects 
can be determined in a natural way only in terms of the 
density matrix. For this reason we transform Eqs. (0) to 
equations for the density matrix <Jij{t) = bi(t)b*{t)- The 

latter is directly related to the amplitudes b{E) via the 
inverse Laplace transform 



+°° dEdE' r 



(27r)2 



h{E)b*iE')e''^^'-^^K 



(8) 
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In order to proceed further with our derivations we have 
to know the i?-dependence of T and A, determined by 
the density of states p{Ea) and the hopping amphtude 
Vt{Ea)^ Eq. (0). We consider below two important cases. 



A. Constant width 

Let us assume that n'^{Ea)p{Ea) is weakly dependent 
on the energy Ea- As a result the width becomes a con- 
stant T{E) = Tq and the energy shift A{E) tends to zero. 
In order to transform Eqs. (0) to the equations for the 
density matrix by using the inverse Laplace transform 
we multiply each of the Eqs. (0) by the correspond- 
ing complex conjugate amplitude h*{E'). For instance by 
multiplying Eq. ( [7a| ) by h^lE') and subtracting the com- 
plex conjugated equation multiplied by bo{E) we obtain 

{E'^E- iTo)b„{E)b*{E') = -t[bo{E) + b*{E')], (9) 

It is quite easy to see that the inverse Laplace transform 
turns this equation to the following one for the den- 
sity matrix 



aoo(i) = -roaoo(t) + Mt) + b*{t)]S{t). 



(10) 



Notice that the initial time in the inverse Laplace trans- 
form (H) corresponds tot = +0, which allows us to ignore 
the term proportional to the 5-function on the right hand 
side of Eq. (0). 

Proceeding in the same way with all the other equa- 
tions , we obtain the following set of equations for the 
density matrix a{t) 



— roa-oo(Oi 

iflaicTaQit) 



(11a) 
(lib) 

(11c) 



with eoa ^ Eq - Ea and aoa = cr*o. 

It is interesting to compare Eqs. (|ll| ) with similar 
Bloch-type equations describing quantum transitions be- 
tween two isolated levels [ p4|Jl(]| ]. In the case of the 
isolated levels {Ei and E2), the equations for the den- 
sity matrix a are symmetric with respect to Ei and i?2. 
Whereas in the case of transition between the isolated 
(Eq) and the continuum of states (Ea) the corresponding 
symmetry, between Eq and Ea , is broken as can be seen, 
for example, in the equation for the off-diagonal term cTao 
where the coupling with aaa is missing. Eq. (11c). 

Solving Eqs. ( pi] ) we find the following expressions for 
the occupation probabilities, ctqo and (Tqq, of the levels 
Eq and E^, respectively pBf : 

CToo(0 =e-^°* , (12a) 

(7aa{t) 



{Ea — EqY ■ 

(^1-2 cos(£'a 



(ro/2)^ 

EQ)te-^"^/'^ + < 



(12b) 



Notice that the line shape, P{Ea) = <Jaa{t — > oo)p, given 
by Eq. (121), is the standard Lorentzian distribution 



P{Ea) 



ro/(27r) 



{Ea — Eq 



(ro/2)^ 



(13) 



with the width Fq corresponding to the inverse life-time 
of the quasi-stationary state, Eq. (12a). 



B. Lorentzian density of states 

Let us consider the same problem studied before but 
with the decay width F in Eq. ( 7a ) dependent on energy. 
This dependence could be generated either by the density 
of states in the reservoir p{Ea) or by hopping amplitude, 
Vta, depending on Ea- Here, we will restrict our selves 
to the frequent situation in which the density of states in 
the reservoir is modulated by a Lorentzian shape, as for 
instance by the resonant cavity. 



P{Ea) 



Fi/27r 



{Ea-E^Y + Tl/A 



(14) 



Here p is a background component and Fi is the width 
of the Lorentzian centered around Ei. Substituting 
Eq. (0) into Eq. ^ and assuming that the hopping am- 
plitude U,{Ea) is weakly dependent on the energy Ea^ 
we can easily evaluate the integral (||). Then Eq. ( |7a| ) 
becomes 



E-Eo+i 



2 E ~ Ei+iTi/2 



bo{E) 



(15) 



with F — 27Til^p . Similar to the previous case we trans- 
form Eqs. (0) into the rate equations for the correspond- 
ing density matrix. For this reason we spilt Eq. ( |l5| ) 
into two coupled equations by introducing an auxiliary 
amplitude 



bi = 



bo. 



(16) 



E- Ei+ iFi/2 
As a result Eqs. (^ become 

E - Eo + i^Ybo{E) - nabi{E) ^ i (17a) 
E - El + i^-ybi{E) - nabo{E) = (17b) 



{E ~ Ea)ba{E) ~ nabo{E) ^ 



(17c) 



The amplitude &i can be interpreted as the probability 
amplitude of finding the electron in the state Ei em- 
bedded in the reservoir. Now using the same procedure 
as that we applied for the derivation of Eqs. (pTl), we 
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transform Eqs. (|T^ into equations for the density ma- 
trix aij (t) : 



O'll(i) = -rifTll + i^ai(Tw - (Toi) 



(18a) 
(18b) 



r + Ti 

o-oi(i) = ieioO-Ql + if^a(croo ~ cril) 7^ o-Ql (18c) 



(18d) 



o'Oq(^) — ^CqoO'oq + «riQ(croo — ctiq) — -^CTQa (18e) 



^lait) = lecifJiQ + iriQ(crio — ctoq) ^^1" 

where eij — Ei — Ej and Cji = Cj*- . 



(18f) 



Although Eqs. (^ are quite different from Eqs.( [ll| ) 
obtained in the previous case, their interpretation is sim- 
ilar. In fact, Eqs. ( p^ represent Bloch-type equations 
describing decay of two-level system, as for instance an 
electron in the coupled double quantum dots, shown in 
Fig. 2. (The latter is described by Eqs. (|l|) for f = 0). 
This is not surprising since the Lorentzian component 
of the density of states corresponds to a resonance state 
embedded in the reservoir. 




Fig. 2: Double-dot system coupled to the reservoir of 
continuum states. The electron is initially in the left 
dot. 



Solving Eqs. (|T|) (or Eqs. (0)) one finds that 
[{E^-E,f+Tl/A] 



oo 



(19) 



Thus, the hne-shape P{Ea) = (Taa{t oo) p{Ea) is not 
a pure Lorentzian one, in contrast with the previous 
case, Eq. (p^. This is a reflection of non-exponential 
decay. The latter is generated by the energy dependence 
of the width T{E) in Eq. (|7a|), as for instance in the 
case of a Lorentzian density of states. Indeed, by solving 
Eqs. (18a)- (18c) for f = we obtain 



CToo(t) 



Ti +UJ + 2i£oi 



2c. -^"P 

Fi - w -I- 2«eoi 
2u~ 



Fi 



■ exp 



, (20) 



where cu = ^y (Fi + 2ieoi)^ — 1617^. Hence, the decay is 
not a pure exponential one. In particular, it follows from 
Eq. ( pO| ) that the probability of finding the electron in 
the initial state for small t is croo(i) = 1 — ^a^^i 



trast with Eq. (12a), describing the exponential decay. 
Notice, that the absence of linear in t term in (Too(i) is a 
prerequisition for the Zeno effect. 

Consider Eq. ( pO| ) for Fi 3> il-a,^oi- One finds that 
the first exponential in Eq. (|2^) dominates the behavior 
of (Too(t) when t increases. Thus 



(Too(t) ^ exp 



4f} 



for t>2/Fi 



(21) 



so that the decay becomes an exponential one. Re- 
spectively, the line-shape P{Ea), Eq. ([l9|), becomes 
Lorentzian in the same limit: 



1 



2tt 



Fi 



(22) 



Using Eqs. (18a)-(18c) one can also evaluate the av- 
eraged decay-time of the electron in the state Eq, given 
by 



nOC pOD 

T = - t croo(i) / croo(t) 

Jo Jo 



(23) 



The last integral can be evaluated directly from 
Eqs. (18a)-(18c). As a result, we obtain for the decay- 
time 



1 + tT 



where 



1 46fo + (f + Fi)^ 
Fi 4r!2(f + Fi) 



(24) 



(25) 



In the case of zero background width, F = 0, the elec- 
tron decay-time is 



T 



1 A{Ei~Eo)^ + Tl 
Fi 4r!2Fi 



(26) 



Hence, instead of the intuitively expected answer 1/Fi, 
we find an enhancement of the decay-time. This enhance- 
ment can be attributed to an oscillatory behavior of the 
electron between the states Eo and Ei, which does not 
exist in the previous case of a constant width. In fact, 
it is even more surprising that for large Fi the electron 
decay-time T increases with Fi. Indeed, in the limit 
Fi ^ Ua,£io the decay-time becomes T ~ Fi/4il2^ in 
accordance with Eqs. (pT]), ([2^). Yet, in the opposite 
limit, eio > Fi, the decay-time is T ~ (1 -f elQ/QD/Ti. 
Thus it is inversely proportional to Fi, as expected. Such 
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an unusual behavior of the decay-time can be understood 
in terms of a broadening of the level Ei due to its cou- 
pling with continuum states Ea, Fig. 2. If Eq = Ei, this 
broadening results in an effective misalignment of the 
levels Eq and Ei , thus destroying the resonant-tunneling 
condition. As a results the decay from the level Eq to con- 
tinuum slows down. On the other hand, if Eq ^ Ei, the 
broadening of the level Ei would effectively diminishes 
the misalignment of the levels i?o,i, so that the decay 
would be accelerated. 

It has been argued that the slowing down of the decay 
rate with Fi in two stage processes, as shown in Fig. 2, 
might be interpreted as the Zeno effect . This would 
imply that the electron decay from the level Ei to the 
states Ea is treated as the measurement of the first tran- 
sition from the level Eq to the level Ei . Such a "measure- 
ment" would localize the electron on the level Eq. Yet, 
the electron in the reservoir is a part of the same mea- 
sured system. Strictly speaking, it cannot be considered 
as an external system to itself. Actually the measure- 
ment must always imply an external macroscopic system 
interacting with the measured electron. Such a measure- 
ment device and its quantum mechanical description is 
provided in the next section. 



III. POINT CONTACT AS A DETECTOR OF AN 
UNSTABLE SYSTEM 
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Fig. 3: The point contact detector near the quantum 
dot. The energy level Eq of the dot is occupied by an 

electron (a), which tunnels from the state i?o, to 
continuous states Ea of the reservoir (b). and /i^j 
are the Fermi levels in the emitter and the collector, 
respectively. 

The dynamics of the entire system is determined 
by the many-body time-dependent Schrodinger equation 
idt\'^{t)) = H\^{t)), where the total Hamiltonian con- 
sists of three components H = Hqd + Hpc + Hint, de- 
scribing the quantum dot, the point contact detector, and 
their mutual interaction, respectively. The first compo- 
nent, Hqd, is given by Eq. (|l|). The two additional com- 
ponents, Hpc and Hint, can also be written in the form 
of a tunneling Hamiltonian as 



As an example of the detector, monitoring decay of an 
unstable system, we consider a point contact placed near 
the quantum dot. Fig. 3. The point contact is connected 
with two separate reservoirs. The reservoirs are filled up 
to the Fermi levels iil and fip, respectively. Therefore 
the current / — e'^TV/2TT flows from the left (emitter) to 
the right reservoir (collector) , where T is the transmission 
coefficient of the point contact and eV = {^l ~ (J-r) is the 
bias voltage (We consider the case of zero temper- 
ature). When the dot is occupied, Fig. 3a, the transmis- 
sion coefficient of the point contact decreases {T' < T) 
due to Coulomb repulsion generated by the electron in 
the dot. Respectively, the current through the quantum 
dot diminishes, I' < I . However, when the electron is 
in the reservoir. Fig. 3b, it is far away from the point 
contact. As a result the transmission coefficient of the 
point contact becomes again T, and the current increases. 
Thus, the point contact does represent a detector, which 
monitors the occupation of the quantum dot. Actually, 
such a point contact detector has been successfully used 
in different experiments |2^. Notice that the current 
variation (/ — /') can be a macroscopic quantity if the 
applied voltage V is large enough. 



Hpc = ^ Eic\ci + ^ Erclcr + ^[fiircjcr + h.c] 

I r l^r 

Hint = '^[Sillrclcr + /i.C.JcJcq , (27) 



where the operators cf{ci) correspond to the creation 
(annihilation) of an electron in state i. The ili^ are the 
hopping amplitudes between the states Ei, E,. of the 
point contact. The quantity S^ir = ^[j. — i^ir repre- 
sents the variation of the point contact hopping ampli- 
tude, when the dot is occupied. 

Consider now the entire system in the initial condi- 
tion corresponding to the occupied quantum dot and 
the reservoirs are filled up to the Fermi levels fiL 
and fiR, Fig. 3a. We denote this state as |5'(0)) = 
Cg|0). This state is not stable: the Hamiltonian (||) re- 
quires it to decay to continuum states having the form 
CiCi • • • |0). In general, the total wave-function 



at the time t can be written as 
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balr{t)cl^clci + 



|0). (28) 



Here b{t) are the probability amplitudes of finding the 
system in the state defined by the corresponding creation 
and annihilation operators. As in the previous case, we 
substitute Eq. ( p8| ) into the time dependent Schrodinger 
equation and use the Laplace transform Eq. (^). Then 
we find an infinite set of algebraic equations for the am- 
plitudes b{E). 



{E - Eo)bo - - ^'ir^ir = I 



l.r 



'Jalr 



(29a) 
(29b) 



{E - Ea)ba - ^abo - ^ illrk 
l,r 

{E + El — Er)bir — ^ ^abalr 

a 

-nlbo^Y.^'i'r~^u'"' =^ (29c) 



I'r' 



These equations look much more complicated than the 
previous ones in Eqs. (0), describing the one-particle de- 
cay. Nevertheless, they can also be transformed into 
Bloch-type rate equations by tracing over the continuous 
degrees of freedom. Such a technique has been derived in 
@yj2|. In this paper we extend this technique by con- 
verting Eqs. (|2^) into the rate equations, even without 
tracing over all continuum states As a results we 

obtain generalized Bloch-type equations that determine 
the energy distribution of the tunneling particles in the 
reservoirs, as well as the time-development of the entire 
system. 

As in the previous case of the one-electron decay we 
replace the sums in Eqs. (|29| ) by the integrals, ~^ 
J p{Ek)dEk. Consider first the terms related to the left 
and the right reservoirs of the detector. Fig. 3. Let us 
assume that the corresponding hopping amplitudes, flir, 
and the density of states p{Ei^r) are weakly dependent on 
the energies: Vlir = V,{Ei,Er) — ft, and p{Ei^r) — Pl.r- 
Let us also assume the large bias limit, eV » ^"^pl^r, so 
the integration limits can be extended to infinity. Then 
the integrations can be treated analytically. Indeed, con- 
sider for instance Eq. ( ^9aD . Replacing the amplitude bir 
by the corresponding expression obtained from Eq. (29c), 
we find 



E-E,- I ^L£l:P^^E!^]bo{E) 



E + El-Er 



(30) 



where J- denotes the other terms where the amplitudes 
bo cannot be factorized out the integral. The integral 



in Eq. (30) is treated in the same way as in Eq. ^ for 
one-electron tunneling. Namely, we split the integral into 
its principle value and singular part. The singular part 
yields iD'/2, where D' = 2nn'^ plPrcV = eT'V/{2TT) 
p9[ , and the principal part is zero. 

Consider now the non-factorizable terms T in Eq. (|3^). 
These terms vanish in the large bias limit. Indeed, all the 
singularities of the amplitude b{E, Ei, Eii, Er, E^i) in the 
£"// -variables lie below the real axis. This can be seen 
directly from Eqs. (|2^) by noting that E lies above the 
real axis in the Laplace transform Eq. (|^). Assuming 
that the transition amplitudes as well as the density 
of states Pl.r are independent of -Bj.r, one can close the 
integration contour in the upper _E;^r-plaiie. Since the 
integrand decreases faster than 1/Ei^r, the resulting in- 
tegral is zero. 

Applying analogous considerations to the other equa- 
tions of the system (|^, we arrive at the following set of 
equations 10|T^ 



E-Eo + i y) ^0 - XI = 

E-E^ + ij^ bc,~nJo = o, 

E + El - Er + i— j blr 



(31a) 
(31b) 



f7'6o - X f^a&a/r - , (31c) 



where D = {2TT)\nir\^ plpr{pi - ^lR) = eTV/{2n). 

Proceeding in the same way with the remaining sum 
over the states Ea we finally arrive to the equations 



D' \ ~ 

E-Eo+i—-I{E)jbo = t, (32a) 

E- E„+i^yb„-nJo = , (32b) 

E + eir + * Y " " ^'^n " ^'^0 ^ ^ ' (^^''^ 



where — Ei — Er and 

TiE) _ f nl(E^)p(E^) 



I{E) = A{E) - 



E~Ea +iD/2 



dEa, (33) 



Let us introduce the reduced density matrix of the en- 
tire system that includes the measured electron and the 
detector: 



a["-"')(i) = (*(i)|n',j)(n,z|vl/(t)) , 



(34) 



where |n, i) denotes the state with n electrons in the right 
reservoir, i.e. J^r'^l'^A^) = '^l"-): and i,j — {0,a} de- 
notes the state of the observed electron. This density 
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matrix can be directly related to the amplitudes b{t) of 

l.r 



Lr 



l.r A 



(35) 



Now we transform Eqs. ( p2| ) into equations for the den- 
sity matrix cTj-" " ^(t), without their explicit solution. It 
can be done by using the same procedure as in Sect. 2, for 
two cases, corresponding to weak and Lorentzian energy 
dependence of the function !(£') , Eq. ( ^3| ) . Since we take 
for the definiteness the amplitude Oq. as independent of 
Ea , these two cases correspond to constant or Lorentzian 
density of states p{Ea) in the reservoir. We demonstrate 
below that the influence of the measurement of the decay 
is distinctly different in both cases. 



A. Constant density of states 

If the product ^'^p{Ea) in Eq. (^) is weakly depen- 
dent on the energy Ea, one can replace r(_E) = Fq, and 
A(£') = in Eqs. (^- Q where Eq = 2T:nlp = const. 
In order to convert Eqs. ( |32| ) into equations for the den- 
sity matrix, we we multiply each of these equations by the 
corresponding complex conjugate amplitude b*{E') and 
subtract the complex co njuga ted equation multiplied by 
ba{E). Eor instance Eq. ( |32b| ) becomes 

{E-E' + iD)b^{E)bl{E') 

= n^ME)b*^{E') - b*{E')bo.{E')] (36) 

Then the inverse Laplace transform, Eq. converts 
Eq. (^6|) to the following equation for the density matrix 

air^ it) = -Dal^^ (t) + [a^ (0 " -^T (01 ■ (37) 

Proceeding in the same way with all other equations 
( p^ ) and integrating over the continuum states of the col- 
lector and the emitter, we obtain the following infinite set 
of equations for the density matrix a{t) 



(") 

00 



-(ro + i?')4"^ 



^00 



D'cr, 



l(n-l) 



00 



v(") 

'aO 



i{Eo 



(") 





D' 



.(") 



(38a) 
(38b) 

(38c) 



(n) _ „(n,n) 
aO 



Here we denoted cr^g' = a 
correspond to o--"'"''(0) 



The initial conditions 



y SioSjoSnoSn'o- Noticc, that 

Eqs. ( p8| ) involves only the diagonal density matrix el- 
ements with respect of the number of electrons in the 



collector n. This, however, does not imply the vanish- 
ing of the off-diagonal terms, cr("'" It means only their 
decoupling from the diagonal terms in the equations of 
motion. This always takes place, in transition between 
continuous states p^ , p8| |. 

Eqs. (^8|) represent a generalization of the previously 
derived Bloch-type equations for quantum transport in 
mesoscopic systems They have clear physical 

interpretation. Consider for instance Eq. (38a) for the 
probability of finding the electron inside the dot and n 
electrons in the collector. This state decays due to one- 
electron hopping to the collector (with the rate D'), or 
due to the electron tunneling out of the dot (with the rate 
Fq). These p roce sses are described by the first ("loss") 
term in Eq. (38a). On the other hand, there exists the 
opposite ("gain") process when the state with (n — 1) 
electrons in the collector converts into the state with n 
electrons in the collector. It also takes place due to pen- 
etration of one electron through the point contact with 
the same rate D'. This process is described by the second 



term in Eq. (|38a|). 

The evolution of the off-diagonal density matrix ele- 



ments cr]^Q {t) is given by Eq. ( p8q ). It can be interpreted 



in the same way as the rate equation for the diagonal 
terms. Notice, however, the difference between the "loss" 
and the "gain" terms. The latter can appear only due 
to coherent transition of the whole linear superposition 

In order to determine the time-evolution of the ob- 
served electron we have to trace out over the detec- 
tor states n in Eqs. (|3^). As a result we obtain the 
following rate equations for the reduced density matrix 

(Too = -Poo-oo (39a) 

&aa — i^a{(^aO — <^0a) (39b) 

ctqo = i{Eo - Ea)(Tao ~ i^aCTQQ a^o ■ (39c) 



Here F^ = {vD — yiy)^ is the decoherence rate, gener- 
ated by the detector Q. 

Let us compare these equations with Eqs. (|ll|), de- 
scribing decay of a single electron, which is not observed 
by the outside detector. We find the di ffere nce o nly i n 
equations for the off-diagonal term, Eqs. (11c) and (39c), 
respectively. The interaction with the detector results in 
an increase of the damping term. All other equations are 
unaffected by the measurement. As a result, the proba- 
bility of finding the system undecayed remains the same 
as in the previous case, croo(i) = exp(— Foi). It means 
that the continuous monitoring of the unstable system 
does not slow down its decay rate. Nevertheless, it af- 
fects the energy distribution of the tunneling electron 
P{Ea) = (Jaa{t — > oo)/9. Indeed, by solving Eqs. (j3^) in 
the limit of i — > oo we find a Lorentzian centered about 
Ea = Eq: 
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(ro + rd)/2^ 



(£;o-So)2 + (ro + rrf)2/4 



(40) 



If there is no coupling with the detector, — 0, the 
Lorentzian width (the hne-width) Tq is the inverse hfe- 



time of the quasi-stationary state, as given by Eq. (13). 
However, as follows from Eq. (pO|), the measurement re- 
sults in a broadening of the line-width due to the de- 
coherence generated by the detector. It now becomes 
To + Td. 

On the first sight this result might look surprising. In- 
deed, it is commonly accepted that the line-width does 
correspond to the life-time. Yet, we demonstrated here 
that it might not be the case, when the system inter- 
acts with an environment (the detector). To understand 
this result, one might think of the following argument. 
Due to the measurement, the energy level Eq suffers an 
additional broadening of the order of r^;. However, this 
broadening does not affect the decay rate of the electron 
Tq, since the exact value of Eq relative to Ea is irrelevant 
to the decay process. In contrast, the probability distri- 
bution P{Ea) is affected because it does depend on the 
position of Eq relative to Ea as it can be seen in Eq. (40). 

Although our result has been proved for a specific de- 
tector, we expect it to be valid for the general case, pro- 
vided that the density of states p and the transition am- 
plitude U,a for the observed electron vary slowly with 
energy. This condition is sufficient to obtain a pure ex- 
ponential decay of the state Eq. On the contrary, if the 
product rj^p(i?c) depends sharply on energy Ea, it yields 
strong £;-dependence of T and A in Eqs. (|^), (|3|). This 
modifies both the exponential dependence of the decay 
probability, croo(i), and the energy distribution, P(Ea), 
as that given by Eqs. ( p9| ) and (pO|). In particular, the 
decay probability for the small t would be proportional to 
t'^. Therefore it is important to investigate the infiucnce 
of the measurement in this case. 



B. Lorentzian density of states 

Consider for the definiteness that the density of states 
p{Ea) of the Lorentzian form, given by Eq. ( p^ ) and the 
amplitude Qa is slowly dependent on the energy Ea- In 
this case one obtains from Eq. d32^ 



^ ' 2 2 E - Ei+i{D + ri)/2 



(41) 



Substituting this result into Eqs. ( |32| ) and introducing 
the auxiliary amplitudes 



bi = 

bllr — 



bo 



E-Ei+i{D + Ti)/2 
E + Ei-Er-Ei+i(D + ri)/2 



Hr 



(42a) 

(42b) 
(42c) 



we can rewrite Eqs. ( p2[ ) in the following form 
.D' + T" 



E-Eo + i 
E- Ei+i 



2 

D + Ti 



bo - Qabi = i , 
bi — ilabo — 



D 



E-Ea+i—j ba-flabo^O , 
D' + f\ ~ 

E + El-Er+ i blr - n'bllr = , 



(43a) 
(43b) 
(43c) 
(43d) 



where f = 27rJ7^p. The amplitudes bi{t), biir{t), . . . are 
the probability amplitudes of finding the electron on the 
resonant level Ei of the continuum, for different numbers 
of electrons in the collector. 

It is quite clear that Eqs. (|43| ) describe an unstable 
two-level system, which is continuously monitored by the 
point contact detector. If F = the entire system is 
equivalent to that shown in Fig. 4, where the penetrabil- 
ity of the point contact is affected only when the observed 
electron occupies the level Eq. 




Fig. 4: The point contact detector near the double-dot. 
The energy level Eq of the first dot is initially occupied 
by an electron. 

Using the same procedure as in the previous case we 
transform Eqs. (^3|) into equations for the density matrix 
(j{t). Finally we obtain the following set of equations, 
which involves only the diagonal terms in n (similar to 
Eqs. (p8|)). For simplicity we consider only the case of 
f = 0. 



■An) 
■'00 



00 



^ '^00 

in) 



Da 



+ iila((Toi' - a 
11 



10 



+ l^a{<Jlo' -<^01 ) 



(44a) 
(44b) 



CToi = ieiQCoi + iila{cr^o - <Jli j 
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D + D' + T 



"01 



(n-l) 



■ (") ■ («) 



D + D' 



.(») 



• (") ■ l"J I -n /' 



(n) 
10 



2D + ri 



(44c) 
(44d) 

(44e) 
(44f) 



These equations describe the microscopic behavior of 
the measured system and the detector at once. Their 
physical meaning is similar to that given for Eqs. (^). In 
order to find the time-evolution of the observed electron 



alone we trace out over the detector states n in Eq. (44). 
As a result, we obtain the following equations describing 
the time evolution of the corresponding reduced density 
matrix elements 



(Too = i^a{croi — CTio) 

till = — FiCTii + ir2a(CTio — CToi) 

coi = leio^oi + i"q(ctoo — ciij (Toi 



O'Oa = ^eaO^OQ + *^^Q(f00 ~ <^lc 



-COo 



(45a) 
(45b) 

(45c) 

(45d) 

(45e) 

(45f) 



where Tfi = {\fD — ^/TyY is the decoherence rate, gen- 
erated by the detector. At first sight Eqs. (^5|) resem- 
ble Eqs. (|l8|), where the background width T is replaced 
by the decoherence width Yd generated by the detector. 
Yet, Yd does not enter the equations for diagonal den- 
sity matrix elements, in cont rast with the background 
width r that appears in Eq. ( 18b ). This very essential 
difference follows from the fact that the measurement is 
a non-invasive one, and therefore it does not result in 
additional dissipation processes. 

Eqs. (^) give a comprehensive description of the mea- 
sured system under continuous monitoring, including its 
energy distribution in the continuum, P{Ea) = (Jaa{t 
oo)p{Ea). The latter can be obtained from Eqs. ( ^ ) 
without their explicit solution. Indeed, let us integrate 
each of Eqs. ( ^5|) in the interval Q <t < oo and take into 
account the initial condition croQ(O) = 1, (7^ (0) = and 
aij{oo) = 0, except for crQQ(oo) ^ 0. Then we obtain the 
following algebraic equations for <J = Jq cr(t)dt: 

iilaio-Qi - CTio) = -1 (46a) 
-riCTii -I- iflaiaio - aoi) = (46b) 

ieioCToi + if^Q(a-oo - o-ii) ctqi = (46c) 

ieaoCToa + iilaiaoo - o-ia) ^^Oa = (46d) 



ieaicria + i^aicria - croa) - -^'^ia = , (46e) 
where aji = a*y The energy distribution is 



IV. ZENO AND ANTI-ZENO EFFECT 



In order to determine how the measurement affects the 
measured system, we solve Eqs. (p5[), ( ^6| ) for ^ 
and compare the solution with the corresponding one, 
obtained for Yd = 0. Let us consider Fi V,a, where 
the system displays an exponential decay, Eq. (|2^), for 
large enough t. 

Let us first examine the probability of finding the sys- 
tem undecayed, (Too(i), given by Eqs. 



f5a|)-(|45c|). This 
quantity as a function of time (in the units fi"^) is shown 
in Fig. 5 in the logarithmic scale for Yi/fia = 10 and (a): 
£10 = ^1-^0 = 0, (b):eio = 10f^„. 



(a) 





10 Slat 



Fig. 5: The probability of finding the electron inside 
the quantum dot at the level Eq (Fig. 4), where (a): 
Eo^ El, and (b): Ei - Eq = 10f7„. The sofid line 
corresponds to F^ = (no measurement) and the 
dashed line to Yd = lOfia- 
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One finds from Fig. 5a that decay rate of tiie elec- 
tron monitored by the detector (the dashed hne) slows 
down, as expected from the Zeno effect. Yet, if the levels 
Eq and El are not aligned, the situation is different. It 
follows from Fig. 5b that the continuous monitoring of 
the decayed system leads to an acceleration of the decay, 
just the opposite to what is expected from Zeno effect 
(the anti-Zeno effect |30 ). However, the later does not 



take place at very short times, where the continuous ob- 
servation still slows down the decay rate. This can be 
seen clearly in Fig. 6, which magnifies the small i-region 
(i < 17-1) of pig^ 5^^^ 



T 



1 46gi + (ri + r,)^ 
Fi mUTi + rd) 



(48) 



This is precisely the inverse exponential factor of Eq. ( |47| ) 
in the limit of Fi 3> i^a- 

Consider now the energy distribution of the tunneling 
electron in the reservoir, P{Ea) = (Tqc(cxd)p(£'q), given 
by Eqs. (^. In the case of ahgned levels, eoi = 0, one 
finds: 



2{ri+Td){TiTd + 4nl)/Tr 



(49) 




Fig. 6: The probability of finding the electron 
undecayed for small t, and {Ei — Eq) = lOfla All 
parameters are the same as in Fig. 5b. 

Actually, the Zeno and anti-Zeno effec ts c a n be re- 
vealed from the analytical solution of Eqs. (|45a| )-(45c) for 
t ^ ■ Indeed, the behavior of croo(t) in this region is 
dominated by the exponent with lesser exponential fac- 
tor. The latter can be fo u nd d irectly from the secular 
determinant of Eqs. (45a)- (45c) (and ctio = (Jqi) in the 
limit of Fi ^ Qa- In which case one obtains 

^oo(0-xp^ <Il±I^.) (47) 



Therefore if eqi < Fi + F^, 



+ (Fi+Fd)2 , 

the second term in the de- 
nominator of Eq. (47) dominates. As a result, the decay 
rate slows down with the decoherence rate F^ (Zeno ef- 
fect). On the other hand, if eoi ^ Fi -t-F^, the first term 
dominates in the denominator of Eq. (^) and the decay 
accelerates with F^; (anti-Zeno effect). 

One can arrive to the same conclusion by evaluating 
the average decay-time T of the electron in the state Eq, 
Eq. (p3|). This can be obtained directly from Eqs. (45a)- 
(45c). In fact we find 



If eoi ^ 0, the analytical expression for P{Ea) is more 
complicated and less transparent. Therefore it is not pre- 
sented here. 

It follows from Eq. ( p9| ) that the width of the energy 
distribution does not correspond anymore to the inverse 
decay time, given by Eqs. (|4^), (p8|). This stays in con- 
trast with Eq. ( p^ ) , where the line- width is precisely the 
inverse decay-time. Moreover the width of the energy 
distribution, given by Eq. (^) always increases with the 
decoherence rate F^, although the corresponding inverse 
decay-time decreases with F^, Eqs. (|47|), (48). Such a 
broadening of the energy distribution is shown explicitly 
in Fig. 7 that displays PiE^), Eq. (||), for F^ = (the 
solid line) and F^ = lOila (the dashed line). Fig. 7 shows 



2 r 



1 . 75 





£a/afi 



Fig. 7: The energy distribution of the tunneling 
electron, P{Ea) for Ei = E'q, and Fi = lOfi^. The solid 
line correspond to F^ = and the dashed line to 
Frf = Fi. 
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that the width of the energy distribution increases very 
strongly with Td, so that P{Ea) is almost flat for Yd ^ 
r^Q. In fact, the same strong broadening of the energy 
distribution takes place for Eq ^ Ei. 

The described above measurement effects can be par- 
tially interpreted in terms of broadening of the level _Eo 
induced by the detector. One can expect that this broad- 
ening would always lead to spreading of the energy dis- 
tribution. On the other hand, its influence on the decay 
rate depends whether the levels Eq and Ei are aligned or 
not. If Eq = El the broadening of the level Eq destroys 
the resonant-tunneling condition, so that the decay to 
continuum slows down. However, if Eq ^ Ei, the same 
broadening would effectively diminish the levels displace- 
ment. As a result, the decay rate should increase. Yet, 
these arguments are not working at very short times, 
when the decay rate slows down even for Eq ^ Ei, as 
shown in Fig. 6. 

In general, such intuitive arguments, based only on the 
level broadening cannot explain all features of the energy 
distribution P{Ea), especially if the coupling of the level 
El with the reservoir is weak. In this case the energy 
distribution, given by Eqs. ( ^ ) shows rather unusual de- 
pendence on the decoherence rate, r^, generated by the 
detector. Let us take, for instance, Eq = 0, £'i = 5ila 
and Fi = CSJIq. The corresponding energy distribu- 
tion P{Ea) is shown in Fig. 8 for F^ = (the solid 
line). Yd — 0.557^ (the dashed line) and F^ — IOJIq. (the 
dashed-dot line). 



1.7 






1 . : 






1.2: 






0.7: 
. 5 


\ 


A 



Fig. 8: The energy distribution P{Ea) as given by 
Eq. (1^) for Eo ^ 0, Ei = 5^^, and Fi = 0.5^2^. The 
solid line correspond to F^ = the dashed line to 
Yd = O.Sr^a, and the dashed-dot line to F^ = lOila- 

If there is no measurement (F^; — 0) the energy dis- 



tribution P{Ea) is strongly peaked near Eq, whereas the 
second peak, at Ea — Ei, is almost invisible. Yet, by 
switching the detector on, the second peak increases with 
Yd- In fact, both peaks are equally pronounced already 
for Fd = Fi. Then, for F^ = 20Fi the second peak be- 
comes a dominant one (the dashed-dot line) , whereas the 
first peak practically disappears. 

Thus we found that continuous monitoring of an elec- 
tron in a double-dot system, weakly coupled with a reser- 
voir changes completely the intensity of spectral lines. 
This can be partially interpreted in the following way. 
Consider first the case of no measurement, Yd = 0. Then 
a small coupling with the reservoir (Fi) cannot essentially 
affect the eigenstates of the system. These are close to 
Eq i, providing that these levels are strongly detuned. 
Therefore the peak in the energy distribution of the tun- 
neling electron lies near Eq, since the electron occupies 
this level. However, if we switch on the detector, it would 
decohere the electron inside the double-dot, before it tun- 
nels to the reservoir. This means that electron tends to 
be equally distributed between the dots. Yet the second 
dot is directly coupled with the reservoir, whereas the 
first dot is not, (Fig. 4). As a result, the energy distri- 
bution would display the peak near Ei. 

V. SUMMARY 

We presented detailed quantum-mechanical analyses 
of a decayed quantum system under continuous monitor- 
ing. As an example we considered an electron tunneling 
from a discrete level of quantum dot to the empty reser- 
voir. The essential point in our analysis is a full quan- 
tum mechanical account of the macroscopic (mesoscopic) 
detector. This allowed us to investigate quantum Zeno 
effect as generated by the Schrodinger evolution of the 
entire system, without invoking the projection postulate. 
In contrast with the previous treatments, which concen- 
trate only on the time-evolution, we analyzed here the 
complementary energy distribution of the decayed sys- 
tem, as well. 

The results of our analysis clearly demonstrate that the 
evolution of a decayed system under continuous monitor- 
ing is crucially related to the energy dependence of the 
density of states in the reservoir and the tunneling am- 
plitude. If these quantities are weakly dependent on the 
energy, a pure exponential decay of the quantum system 
takes place. In this case the measurement does not affect 
the decay rate. Yet, the energy distribution of the elec- 
tron in the reservoir becomes strongly broadened. As a 
result, the corresponding line-width is no more given by 
the inverse decay time. 

In the case of strong energy dependence of the density 
of states (or of the tunneling amplitude) the situation is 
different. As an example we considered the Lorentzian 
component in the density of states of the width greater 
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than the tunnchng amphtudc. In this case the decay is 
not an exponential one for small times. But it becomes 
exponential when the time increases. Nevertheless, in 
contrast with the previous case the measurement affects 
the decay rate, even in the exponential regime. The ef- 
fect, however, depends on displacement between the ini- 
tial state energy and the mid-energy of the Lorcntzian 
density of states. If these energies are close to each other, 
we obtain the expected Zeno effect, i.e. the slowing down 
of the decay rate at any time. Yet, in the opposite case, 
when the displacement of these energies is larger than the 
Lorentzian width, the measurement induces an acceler- 
ation of the decay rate in the regime of the exponential 
decay (the anti-Zeno effect). At small times, however, 
the measurement always slows down the decay rate. 

With respect to the energy distribution of the tunnel- 
ing particle, the measurement always generates strong 
broadening. However, for very narrow Lorentzian shape 
of the density of states (the Lorentzian width is smaller 
than the tunneling amplitude) the effect of the measure- 
ment is more peculiar. This takes place where the dis- 
placement between the initial energy and the Lorentzian 
mid-energy is large. Then, if there is no measurement, 
the energy distribution peaks near the initial state en- 
ergy, whereas the second peak (near the Lorentzian cen- 
ter) is very weak. However, when the measurement is 
switched on, the peak at the initial energy disappears 
with the increase of decoherence, generated by the de- 
tector. On the other hand, the second peak near the 
Lorentzian mid-energy arises. 

All these results were obtained without taking into ac- 
count the coupling of the fermionic quantum dots to the 
radiation or to the phonon fields. In fact, in a weak 
coupling regime our treatment can be extended to this 
case. Yet, we do not expect any essential modification of 
the results obtained in this paper. The strong coupling 
regime, however, needs a special consideration. 

Our analysis shows different effects generated by con- 
tinuous measurement. All of them were obtained solely 
from the time-dependent Schrodinger equation applied 
to the entire system. Moreover, it looks that some of 
these effects might be in a contradiction with the projec- 
tion postulate, which leads to the slowing of the decay 
rate. Yet, this this point deserves special investigation 
on the level of the macroscopic description of the mea- 
suring device. It would involve a proper definition of the 
measurement time, needed for application of the projec- 
tion postulate. Although we are not coming to this point 
in our paper, we nevertheless consider our analysis as a 
necessary step for a better understanding of the measure- 
ment problem and the nature of the projection postulate. 
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